Towards a wave-extraction method for numerical relativity: III. Analytical examples 

for the Beetle-Burko radiation scalar 



Lior M. Burko 1,2 , Thomas W. Baumgarte 3,4 , and Christopher Beetle 5 ' 6 
1 Department of Physics and Astronomy, Bates College, Lewiston, Maine 0424-0 
2 Department of Physics, University of Alabama in Huntsville, Huntsville, Alabama 35899 
3 Department of Physics and Astronomy, Bowdoin College, Brunswick, Maine 04011 
4 Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801 
Department of Physics, Florida Atlantic University, Boca Raton, Florida 33431 
6 Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada 
(Dated: draft of February 7, 2008) 

Beetle and Burko recently introduced a background-independent scalar curvature invariant for 
general relativity that carries information about the gravitational radiation in generic spacetimes, in 
cases where such radiation is incontrovertibly defined. In this paper we adopt a formalism that only 
uses spatial data as they are used in numerical relativity and compute the Beetle-Burko radiation 
scalar for a number of analytical examples, specifically linearized Einstein-Rosen cylindrical waves, 
linearized quadrupole waves, the Kerr spacetime, Bowen-York initial data, and the Kasner space- 
time. These examples illustrate how the Beetle-Burko radiation scalar can be used to examine the 
gravitational wave content of numerically generated spacetimes, and how it may provide a useful 
diagnostic for initial data sets. 

PACS numbers: 04.25.Dm, 04.30.Nk, 04.70.Bw 

I. INTRODUCTION 

Compact binary systems of co-orbiting neutron stars or black holes are among the most promising sources of 
gravitational waves for a set of recently-constructed interferometric detectors, including the Laser Interferometer 
Gravitational Wave Observatory (LIGO). A substantial effort has been underway for some years now to model 
these systems theoretically, and thereby establish detailed predictions for comparison with experimental data. These 
predictions will likely be needed not only for verification of the theory, but for constructing filters needed to separate 
signal from noise in the detectors. Because a compact binary system admits no obvious symmetry or approximation, 
the theoretical problem is being tackled primarily within the field of numerical relativity, which aims to integrate 
directly the full Einstein field equations on the computer. 

A numerical solution of the Einstein equations usually proceeds in three steps (see, e.g., [l| for a recent review and 
references) . The first step solves the constraint equations of general relativity to find initial data for the problem. The 
second then evolves these data forward in time using the remaining field equations. In the third step the numerical 
solution of the second step is interpreted physically. This paper — as the other papers in this series — is interested in 
the third, last step of numerical relativity. Although conceptually straightforward, this program encounters many 
technical difficulties which make it seem unlikely that numerical relativity will be able to evolve compact binaries for 
many orbits (see, e.g., [HE1EI and p| for recent progress). Within the field, the consensus answer to these difficulties 
is to start the simulations using initial data which describe the system at small binary separations, shortly before 
coalescence. This approach introduces a fresh conceptual problem, however. It is far from clear how to pick initial 
data at small separation so as to replicate the outgoing wave pattern which would arise during a gradual inspiral 
from large separation. Intuitively, generic data describing a close binary will contain spurious radiation originating in 
the particular mathematical technique used to solve the constraint equations rather than in the astrophysics of the 
problem. In fact, it is known that different methods of solving the constraint equations lead to physically distinct 
initial data (see, e.g., [UlaO])- 

This paper is the third in a series which outlines a scheme to analyze the radiation content of a given spacetime 
using invariant techniques based solely on the physical metric. More specifically, this paper uses an invariant radiation 
scalar to measure the amount of radiation, including spurious radiation, in a variety of initial data sets familiar to 
the numerical relativity community. In general, a radiation scalar is a scalar function on spacetime derived from 
the physical metric which, when spacetime unambiguously describes gravitational radiation propagating in some 
background, manifestly contains information only about the radiative degrees of freedom. There may be many such 
metric invariants which could be defined. We will focus on just one, denoted £, which we call the BB radiation 
scalar |8j- This particular scalar seems well-suited to problems, like binary black- hole inspiral, where the quiescent 
end-state describes weak radiation propagating outside a single black hole approaching equilibrium. Ultimately, the 
idea would be to minimize in some sense this measure of radiation content and thereby identify more physically 
reasonable binary initial data sets. While this intuitive idea is easy to state, there are a number of caveats to taking 
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the picture it suggests too seriously. Nonetheless, the results of this paper suggest the radiation scalar approach could 
used effectively in this way. 

The definition of the BB scalar, and the scheme on which it is based, has certain strengths and limitations. In this 
context, the key strength of this approach lies in its direct applicability to initial data sets. (In addition, the BB scalar 
may also have intriguing applications for evolution data.) The only currently-known unambiguous way to minimize 
the spurious radiation content of such data is to use the evolution equations themselves. As evolution proceeds, 
spurious radiation will naturally propagate outward, and eventually off the numerical domain for the problem. The 
gravitational radiation remaining in the instantaneous data at later times would then be the astrophysically sound 
radiation generated by the binary at earlier times. However, the dynamical evolution also introduces numerical error, 
and the difficulties associated with this dynamical evolution, especially for binary black hole systems, motivated the 
construction of small-separation initial data in the first place. In contrast, the BB scalar may be calculated directly 
from an initial data set (see the first paper in this series hereafter Paper I). No evolution is required. The approach 
also offers a potentially useful tool for the extraction of gravitational radiation from numerically evolved spacetimes 
(see hereafter Paper II). However, here we focus exclusively on the application to initial data sets. 

The main risk of the BB scalar approach lies in trying to push its physical interpretation beyond its appropriate 
scope. This point has been emphasized in previous papers 8,., 9, 10] in this series. While mainly conceptual in origin, 
this issue raises genuine practical problems in applications like those considered in this paper. 

The interpretation of the BB scalar is unambiguous when spacetime describes weak-wave perturbations of a single 
black hole, an approximation which becomes increasingly accurate at late times after coalescence of a binary system. 
In particular, it is certain to be small in this perturbative context. However, it remains calculable much more generally, 
although this precise physical interpretation is lost. This is hardly surprising since even the notion of gravitational 
radiation cannot be defined unambiguously in general regions of spacetime. More pragmatically, there is nothing to 
guarantee that the radiation scalar cannot be small in this broader context, where spacetime might intuitively be 
said to describe strong radiation fields. That is, although the radiation scalar is definitely small when little radiation 
is present, it may also be small coincident ally even when this is not the case. This can happen because, whereas 
ingoing and outgoing gravitational radiation in perturbation theory is intuitively associated with the two Newman- 
Penrose Weyl curvature components tpo an d tpi measured in a particular null-tetrad frame, the BB scalar measures 
just their product: £ = ipoipi- Thus, in an extreme example, even though the outgoing component of the radiation 
field described by i/>4 may be non-perturbatively large, if the ingoing component ipo vanishes then so does £. However, 
this cannot happen accidentally. Intuitively, outgoing radiation will scatter off the fixed potential described by the 
background in which it propagates to produce an ingoing signal. The outgoing signal would have to be very carefully 
tuned to the background in order to have the the net effects of this back-scattering cancel exactly throughout some 
region of spacetime. The mathematical expression of this fact lies in the Tcukolsky-Starobinsky identities, which 
interrelate the derivatives of the two seemingly independent components of the radiation field ipQ and ip^ in a Kerr 
black-hole background. Following this extreme example, one might reasonably expect situations wherein the BB 
scalar is coincidentally small are rare and non-generic. This paper tests that expectation using exact, analytical 
initial data sets whose gravitational wave content is already — independently — rather well understood. The results are 
encouraging. We find that, at least in these specific cases, the expectation that the radiation scalar should not be 
small coincidentally in typical situations of practical interest is indeed supported by the facts. 

This paper is organized as follows. In Section [HJ we summarize the recipe developed in Paper I for extracting 
the BB scalar from 3+1 Cauchy data for the gravitational field. The following three sections explicitly calculate 
the BB scalar in three different continuous families of analytical initial data sets. Section ITTT1 focuses on linearized 
wave propagating on a flat background. This is a worthwhile test because the BB scalar was explicitly developed as 
a tool to analyze such linearized waves propagating around a Kerr black hole. In fact, the presence of a dominant 
background component, associated with the "Coulombic" Kerr geometry, has been used explicitly to define the BB 
scalar. It is therefore noteworthy that the approach is also effective when such a dominant Coulombic field is absent. 
Section llVl turns to genuine black hole solutions, analyzing Bowen-York initial data for single black holes using the 
BB scalar. This is probably the most physically significant application considered here. It finds that, whereas the BB 
scalar does vanish in the Kerr spacetime where clearly no radiation is present, it does not vanish for the Bowen-York 
data. We explicitly find a quartic growth of £ with the angular momentum L of the Bowen-York hole. This is exactly 
what one would naively have expected. Section calculates the BB scalar in the cosmological Kasner solution. We 
do this to illustrate both the limitations and the ancillary benefits of our approach. The BB scalar does not vanish in 
these spacetimes, even though they are understood to contain no radiation. However, this is not surprising since these 
cosmological solutions do not contain a radiation zone in which the usual physical meaning could be attached to the 
BB scalar. Nonetheless, we find that it can be used as an interesting gauge-invariant tool to partially characterize the 
gravitational field in this context. One must only resist the temptation to think of it as a radiation scalar. Finally, 
Section IVII summarizes our results and makes a few further comments regarding their interpretation. 
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II. RECIPE 

This Section closely follows the method of paper I to construct the BB scalar £ from numerical relativity data. We 
refer to paper I [t| for further details and discussion of this recipe. It is recapitulated here only for convenience. 

We begin our discussion at the level of spacetime geometry, where the Weyl curvature tensor picks out exactly 
three null tetrads, up to certain scaling transformations, at a generic event 1 . These are the "transverse frames" 
in which the Newman-Penrose curvature components ipi and ip3 vanish. At sufficiently large distances from the 
sources, in a radiation zone, it is possible to single out one of these three, called the quasi-Kinnersley frame, on 
which our definitions are based. The name derives from a particular null tetrad introduced by Kinnersley on a 
particular spacetime, the stationary Kerr black hole, and subsequently used by Teukolsky to analyze the gravitational 
radiation content of perturbed black hole spacetimes. Teukolsky found that, when one makes perturbative corrections 
to the background Kinnersley tetrad to maintain vanishing Newman-Penrose curvature components ipi and 03 in the 
perturbed spacetime, the gravitational radiation is entirely encoded in the two curvature components ipo an d 04 (see 
also H3)- The remaining non-zero component, ifo, is intuitively associated with a Coulombic part of the field, and is 
dominated by terms deriving from the stationary background. 

Once the quasi-Kinnersley frame has been identified in a radiation zone, its definition may then be carried into 
more general regions of spacetime using a continuity principle. The quasi-Kinnersley frame is always one of the three 
transverse frames at any point of spacetime, whether that point lies in a radiation zone or not. In these transverse 
frames, the Newman-Penrose scalar 02 will typically take three distinct values, only one of which will lead to a 
continuous function when compared with known quasi-Kinnersley values of ip2 at nearby points. The frame so chosen 
is the quasi-Kinnersley frame. Put another way, the key idea is that there is no ambiguity in identifying the quasi- 
Kinnersley frame in a radiation zone, and maintaining continuity everywhere in spacetime generally determines a 
unique "transverse frame field" throughout. This is the quasi-Kinnersley frame, defined globally. The results of this 
paper will indicate that there ought to be little doubt in practice concerning which value of ip2 will give the desired 
continuous extension. See the Discussion for further comments. 

This paper will consider two curvature scalars defined using the quasi-Kinnersley frame: 

X := 02 and £ := ip 4 . (1) 

Each of these functions is understood to be defined on spacetime by evaluating the indicated Newman-Penrose 
curvature components in the quasi-Kinnersley frame. The first of these quantities, x, is called the Coulomb scalar. 
The nomenclature derives from intuition in the Teukolsky formalism: -02 is dominated by the stationary, Coulombic 
field of the background Kerr black hole. The second quantity, £, is the BB radiation scalar. In a radiation zone 
accurately described by the Teukolsky formalism, it manifestly depends only on the radiative degrees of freedom. 
Although the components 0o and 04 are not separately preserved by the residual scaling transformations of the null 
tetrad allowed by the Weyl curvature, their product is. Thus, the BB scalar is a true invariant derived from spacetime 
curvature. 

Although certainly not a radiation scalar, the Coulomb scalar is still quite useful in our approach. First, as 
described above, demanding continuity of x plays the key role in defining the quasi-Kinnersley frame. We could also 
have used continuity of £ to do this, but £ is quadratic in the curvature while whereas x is only linear. The radiation 
scalar therefore falls of much more rapidly (£ ~ r -6 ) at large distance than the Coulomb scalar (x ~ r~ 3 ), making 
the numerical problem of determining continuity much harder. Second, the Coulomb scalar can be used to extract 
information about the single black hole underlying a perturbation-theoretic description of a given spacetime. Thus, 
for example, it may be used at late times to determine, or at least constrain estimates of, the mass and spin of the 
final black hole resulting from coalescence in a background-independent way. 

Numerical relativity is not primarily concerned with spacetime geometry. Rather, it solves a Cauchy problem for 
general relativity, calculating a series of spatial geometries parameterized by a fiducial time variable. These may then 
be assembled to form spacetime. The assembly procedure can of course be done analytically, and used to express the 
BB scalar in terms of given Cauchy data. This has been done in Paper I, and we describe the results here. 

In a 3 + 1 decomposition, the gravitational fields are expressed in terms of a spatial metric 7^ and an extrinsic 
curvature Kij. The electric and magnetic components Eij and Bij of the spacetime Weyl tensor can be computed 



Where spacetime is of algebraic Petrov types I or D, there are three such frames. Where spacetime is type II, there is only one. In 
the remaining cases, types III and N, no such frames exist. These last cases, though, are far from generic, and ought not arise in the 
situations of interest here. By far, the most probable algebraic type at a typical point of the spacetimes we study is type I. 
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from these data using 

Eij = II,, • A A, , K.i K; \-S, s (2) 

Eij = £ij - ^ij^Su (3) 

where Sij = 7j°Sfj Tab is the spatial projection of the stress energy tensor T a f,, and i?y is the spatial Ricci tensor 
(see, e.g., Ref. |l| for more detail). The definition for the magnetic part of the spacetime Weyl tensor depends on the 
convention used for the extrinsic curvature and for the spacetime volume element. We adopt here the convention that 



K ij = ~ \tn 7ij . (4) 

where C n denotes the Lie derivative in the direction of the normal n a to the 3-hypersurface. Our convention for the 
spacetime volume element is 

£abcd = 24i£ [a n 6 m c m d ] (5) 

(see paper I for more detail.) We next define the magnetic part of the spacetime Weyl tensor 

Btf = -e^VfcAy (6) 
Bij = Bfa). (7) 

By construction, Eij and are both symmetric and traceless. We then define the complex tensor 2 

C i j =E i j +iB i j , (8) 

from which the scalar curvature invariants / and J can be computed as 



I = ^ (9) 



and 

1 
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J = --C i j C j l C\. (10) 



Once a basis is chosen, the tensor C % j becomes simply a 3 x 3 matrix. In ordinary linear algebra, it is well-known 
that any scalar invariant which may be extracted from such a matrix must be expressible as a function of three basic 
invariants. These basic invariants may be chosen to be the eigenvalues or, for example, the traces of the first three 
powers of the matrix in question. The trace of our matrix vanishes because it derives from the spacetime Weyl tensor, 
and the traces of its square and cube are basically / and J. Thus, we should expect our Coulomb and BB scalars may 
be expressed as functions of I and J. To do this, we recall the Baker-Campanelli speciality index fl3j defined by 

S:=27^. (11) 

In terms of S 1 , the Coulomb scalar can be written as 

o,±_ M W^W + WxiS)- 1 ' 3 



2 Note that the relative sign in the definition of C 1 j here differs from that used in Paper I. At first glance, it would appear that the 
quantities we calculate here are therefore the complex conjugates of those calculated there. The difference, however, is illusory. Our 
previous papers have tacitly assumed a definition for the extrinsic curvature with the opposite sign to that in Eq. J1J. While the previous 
convention is popular in the mathematical relativity community (see, e.g., Il 21 ). the current one is used more often in the numerical 
relativity community. We therefore adopt it here for the convenience of the latter. As a result of this shift of convention, the sign of the 
tensor Bij defined above is reversed. We compensate this reversal in the definition of C l j so that this tensor takes the exact same value 
on a Cauchy slice with our new conventions as it did with the old. With this one modification, all other formulae from the previous 
papers continue to hold unmodified. 
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where W X (S) = </S - VS - 1, while the BB scalar i 



is 



* 4 



W^(S) 1/3 ~w^s)- 1/3 



(13) 



where W^(S) = 2S — 1 + 2y/S(S — 1). In general, both the Coulomb and BB scalars admit three distinct complex 
roots, associated with the three distinct transverse frames. We denote these three roots with the superscripts "0,±". 
The value in the principal branch of the root functions above has superscript 0, and this branch defines the quasi- 
Kinnersley frame in an appropriate radiation zone. The other two branches have superscripts ±. Either of these may 
be labeled the quasi-Kinncrslcy value in the interior of spacetime, depending on what is needed to maintain continuity 
while coming in from the radiation zone. We will illustrate this point explicitly below. 

In addition to the direct calculation above, there are two other ways to compute £ 0,± and x 0,± - The first is to solve 
the eigenvalue problem for the complex spatial tensor C*. : 



o,±\ 



(14) 



As shown in Paper I, the three eigenvalues A turn out to be exactly twice the Weyl scalar fa evaluated in the three 
transverse frames. Thus, the Coulomb scalar at any point of spacetime is always half of one of the eigenvalues of C^. 
If these eigenvalues are calculated, this approach does not say which is which, so it is not immediately clear which 
eigenvalue is associated with the principal branch of i|12fl , and therefore with the quasi-Kinnersley frame. However, it 
has been shown in Paper II that, among the three branches, the principal value has the largest modulus in a radiation 
zone. Thus, the prescription to calculate the Coulomb scalar is to take half the eigenvalue of largest modulus at large 
distances, and extend inward to strong-field regions always choosing the unique eigenvalue which keeps the function 
X smooth. In a transverse frame, the Newman-Penrose Weyl scalars fa and fa vanish by definition, and this gives a 
particularly simple relation between the Coulomb and BB scalars: 



o.± 



(fafa) ' 



(15) 



Thus, if the Coulomb scalar is known, the radiation scalar follows immediately. The second way to calculate the 
Coulomb and BB scalars is to use series expansions of the exact formulae above. As shown in paper I, for a spacetime 
that asymptotically is perturbed about algebraic speciality, we can expand 



X 



J 



1 



and 



£ = - 5 (S-D 



4 (S _ 1) + |L (S _ 1)2 + 



l--(5-l) + — (5-1) 2 . 
27 v ' 729 v ' 



(16) 



(17) 



These expansions explicitly use the principal values of the root functions in (|12(l and (jl3|l . and therefore allow us to 
identify the quasi-Kinnersley \ an d £ asymptotically . This technique can be useful numerically when all three values 



,o,± 



become small. We illustrate this fact in Sections IIV Bl for Bowen-York data and IVl for the Kasner spacetime. 



III. LINEARIZED WAVES ON A FLAT BACKGROUND 



Probably the first example that comes to mind when evaluating a radiation scalar are pure wave solutions. Unfor- 
tunately, the simplest such solution, plane waves, is too simple to be of interest. Plane waves are Petrov type N, and 
hence algebraically special. In a type N spacetime, a frame can always be found so that fa is the only non-vanishing 
Weyl scalar. This frame is automatically transverse, so that £ = fafa vanishes identically. As the BB scalar is frame- 
independent, it vanishes also in all other, non-canonical frames. In a plane wave spacetime the BB scalar vanishes 
despite the presence of gravitational radiation. 

We therefore turn to two less trivial examples, namely linearized time-symmetric Einstein-Rosen cylindrical gravi- 
tational waves in Section lill Al and linearized even-parity quadrupole Teukolsky waves in Section Till Bl 



A. Linearized Einstein— Rosen cylindrical waves 



The metric for Einstein-Rosen cylindrical waves of amplitude B and pulse width a |l4l ] is given by 

dg 2 = e 2[ 7 ( t ,p)-<A( M) ] (_ dt 2 + dp 2^ + e ^p) dz l + p 2 e -2W, P ) d tf t 



(18) 
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where 



^(a + it) 2 +p 2 y/(a - it) 2 + p 2 



(19) 



(which is a real function.) In what follows we consider linearized waves, so that we can take the function j(t,p), 
which is quadratic in B, to vanish. The solution is time-symmetric about t — 0, and describes cylindrical waves that 
are imploding on the symmetry axis for t < and exploding for t > 0. 

In the context of numerical relativity, the spacetime metric is usually expressed in terms of a lapse a, a shift (3 l and 
a spatial metric jij . These quantities can be found be identifying the spacetime metric IjlSf) with the "ADM" form 

ds 2 = -a 2 dt 2 + ^j(dx* + f3 l dt){dx 3 + ftdt) , (20) 

Here the shift vanishes identically, (3 l — 0, the lapse is 

a = exp[-i/>(t,p)], (21) 

and the spatial metric is 

d(J 2 = e 2 HKt,,)] dp 2 + e W(t,P) dz 2 + p 2 e -W(t,P) d( p 2 . (22) 
We also compute the extrinsic curvature Kij from Eq. Q 

Kij = - ^ (d^ - DiPj - DjPi) , (23) 

where Di is the covariant derivative compatible with the spatial metric 7^, and where we use the sign convention 
typically used in numerical relativity. 

At the moment of time symmetry t = 0, K%j vanishes identically, which greatly simplifies the problem. The only 
non- vanishing term in the electric © or magnetic Q components of the Weyl tensor is then the Ricci tensor Rij in 
the electric part Eij , which yields 

E" p = -2B {p 2 + a 2 Y m 

E% = 2B (2a 2 -p 2 ) (p* + a 2 y 5/2 

E\ = 2B (2p 2 -a 2 ) (p 2 + a 2 y b/2 (24) 

to linear order in B (all other components vanish.) At t = we therefore have C^- = E z j. Using JjJJ and (|10|) we 
find the curvature invariants 

T 1 o P 2 q4 Z a2 P 2 + P 4 f>>x\ 

1 = UB (P 2 + a 2 f (25) 



and 



and from (|llfl the specialty index 



(2p 2 - a 2 ) 2 (p 2 -2a 2 ) 2 (p 2 + a 2 ) 2 

4(a 4 - a 2 p 2 + p 4 ) 3 ' 1 j 

We plot the specialty index S as a function of p in the upper panel of Fig. Interestingly, S is independent of the 
wave amplitude B and therefore does not provide a measure of the strength of the gravitational radiation. 

Since S is independent of B, the Coulomb scalar \ 1121) must scale with B and the BB scalar £ (|13fl with i3 2 . 
Reversing this argument we can see that the speciality index S must be independent of the wave amplitude whenever 
the Coulomb scalar is perturbative and scales with B, while the BB scalar scales with B 2 . We find the same behavior 
for the Teukolsky waves in Section lill Bl 

For the remainder of this Section we determine the Coulomb scalar \ and the BB radiation scalar £ explicitly as 
functions of p, adopting the following strategy. There are three regions for which the spacetime is close to speciality, 
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i.e., for which S — 1 -C 1: p S> a, p — a and p <C a (compare Fig. QJ) We then identify x with an eigenvalue of C l j. 
Specifically, we require that at asymptotically great distances (p a) the expansion holds, and then identify x 
with an eigenvalue of C l - — which is equivalent to picking a particular branch in the cubic roots in (|12|1 . From \ we 
can then compute £ using (|15|l . Finally, we require continuity and differentiability as we move to smaller and smaller 
values of p/a. 

Before we proceed we list the three eigenvalues of Cj, which we denote A p , A 2 and A^, 

A p = -2B (p 2 +a 2 y 3/2 

\ z - 2B (2a 2 -p 1 ) (p 2 +a 2 y 5/2 

A = 2B (2p 2 -a 2 ) (p 2 +a 2 y 5/2 . (28) 

Notice that I = (1/2) £\ \ 2 and J = (l/2)n,A l as expected. 

To find x and £ we examine the asymptotic region p ^> a. In this limit, 

I-^rB 2 p»a, (29) 

and 

J ~ -4r^ 3 p> a, (30) 

so that (fT?))l yields 

X~-3y~^B. <°» a ' ( 31 ) 

This the quasi-Kinnersley x- We compare this result with the limiting values of the three eigenvalues (|28|) of the 
tensor C*., 

A p PS —-^t -B, A z m — ^r-B, A^^B p>a, (32) 

pd 

The eigenvalue that corresponds to x is the one that satisfies A = 2%, and therefore we conclude that for p ^> a, the 
correct eigenvalue is A^. (Alternatively, this is the eigenvalue with the greatest modulus — see paper II.) Next, we 
require that x is smooth to extend x to other regions of the spacetime. In principle, we could switch to a different 
eigenvalue at a branch point. This would make x discontinuous, however, which is not desirable. A distributional x 
would suggest similar characteristics of the physical Weyl scalar, which is clearly not the case in this example. The 
only identification for £ that makes it smooth throughout is that x = A^/2 not just for p a, but everywhere. We 
therefore have 

X= WTaW B (33) 
everywhere. We can now compute the BB radiation scalar £ from (|15|l . which yields 

e = I-3x 2 = 9(^)5^- (34) 

Our results for x an d £ are plotted in Fig. \T]m the middle and lower panels, correspondingly. 

In Appendix^]we present an alternative, yet equivalent, derivation of x an d £, that is based directly on Eqs. I|12|) 
and l(T3|l . respectively. 



B. Linearized quadrupole waves 



Einstein Rosen cylindrical waves are not asymptotically flat, which limits their usefulness as testbeds for numerical 
relativity. (In addition, Einstein-Rosen cylindrical waves do not represent the exterior solution of bounded radiating 
systems.) We therefore briefly discuss a second example of linearized gravitational waves over a flat background, 
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FIG. 1: The speciality index 5* (upper panel), the dimensionless Coulomb scalar a 3 \/B (middle panel), and the dimensionless 
BB scalar a 6 £/ B 2 (lower panel), as functions of p/a for time-symmetric linearized Einstein-Rosen cylindrical waves of pulse 
width a, at t = 0. 



namely Teukolsky's even-parity linearized quadrupole gravitational waves in vacuum with azimuthal symmetry |15| . 
This spacetime is algebraically-general (Petrov type I), that asymptotically at great distances (A 1 / 2 r>l) approaches 
algebraic speciality (Petrov type III). This algebraic class will be important for the interpretation of the speciality 
index S and the BB scalar £. 

From the spacetime metric (Eq. (5) in p^| ) we find the lapse a to be unity, the shift j3 l to vanish, and we can 
identify the spatial metric 7^ as 

da 2 = [1 + (2 -3 sin 2 6)A]dr 2 ~6Br sin 2 9 cos 6 dr dO + {I - A + 3 sin 2 6C)r 2 d9 2 

+ [1 + (3sin 2 6»- l)A- 3 sin 2 (9C]r 2 sin 2 6>d</> 2 . (35) 

Here the coefficients A, B and C are given in terms of a function F(x) := F(t =p r) and its derivatives F^ n ' = 
d n F{x)/ dx n \ x=tTri 

/ F m 3F (i) 3F \ 

- 3 — ± — + _ (36) 



F (3) 3i?(2) 6F (1) Q F ^ 

®t = -(±— + — ± — + ^) (37) 



1 /FW 2F<® 9F( 2 ) 21F« 21F\ , . 

C t = 5 (— ±^ + ^±^^ + — J- (38) 

where the upper/lower signs correspond to that in the definition of F. We would like to construct time-symmetric 
data with Kij — at t = 0, and therefore choose linear combinations A = A- — A+, B = B- — B+, and C = C_ — C+. 
For concreteness, we choose the function 

F(x) = e{t±rf e- xi - t±r ^ 2 (39) 
where e is a measure of the wave amplitude and A is related to the wavelength. 
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As for Einstein-Rosen waves, the vanishing of the extrinsic curvature at t = greatly simplifies the problem. Since 
the spatial Ricci tensor R^ is linear in e, we have S^j = R^ to leading order in e and also fiy = 0. This scaling is 
sufficient to determine the scaling of the BB radiation scalar £: Since = E l j is linear in e, I must scale with e 2 
and J with e 3 . We again find that the speciality index S is independent of e. ft does depend on the choice of F(x), 
but is independent of the wave amplitude. That is, S is insensitive to the amplitude of the waves. As we demonstrate 
in Appendix El the speciality index S — > as A 1 / 2 r — > oo. However, as spacetime is asymptotically neither Petrov 
types nor N, while I, J — > 0, it is algebraically special and of Petrov type III. This property is not accounted for by 
the speciality index S, which does not approach unity as A 1 / 2 r — * oo. 

The Coulomb scalar \ again scales with the wave amplitude parameter e, while the BB scalar £ scales with its 
square, e 2 . 

Explicit expressions for the Coulomb and BB radiation scalar can be constructed with an analysis similar to that 
in Section IIII Al for Einstein-Rosen waves. The calculations are significantly more involved, however, and do not 
necessarily provide new insight. Instead of going through these calculations we assure the reader that the results are 
qualitatively very similar to those for Einstein-Rosen waves, and refer to Appendix [B] for some details. 



IV. ROTATING BLACK HOLES 



A. Kerr black holes 



The Kerr solution is of Petrov type D and hence algebraically special with 

5=1, (40) 

independently of the black hole angular momentum L = aM. From Eq. H17|) it is evident that the BB radiation scalar 
vanishes, £ = 0. This can also be seen by observing that Wj = 1 when S = 1. The three cubic roots of unity are 
l,e 2l7r ' 3 and e 4l,r / 3 , which, when inserted into £ 0,± , yield zero and twice 37/4. These are the three roots found in 
paper II, and the BB radiation scalar corresponds to the vanishing value. Given that the Kerr spacetime is stationary 
and does not contain any gravitational radiation it is reassuring that the radiation scalar vanishes. 

Even though we already know the result, it is useful to demonstrate that the recipe laid out in Section[n]would lead 
to the same conclusion. To do so we start with the spacetime metric of a Kerr black hole, which in Boyer-Lindquist 
coordinates can be written as 

£ p 2 



ds z = -- dt 2 + ^sm z 9 \d(j> dt) + '— dR z + p z dQ* (41) 



where we have used the abbreviations A = R 2 — 2MR + a 2 , p 2 = R 2 + a 2 cos 2 9, and £ = (R 2 + a 2 )(R 2 + a 2 cos 
2a 2 MR sin 2 0. As in Section UlI Al we can identify the spatial metric of constant Boyer-Linquist time slices as 



n 2 S 

~iijdx x dx? = '--dR 2 + p 2 d9 2 + ^ sin 2 6d(f> 2 , (42) 
A p z 



the lapse as 



and the shift vector as 



* - (^) 



(3 a = (0, 0, -2MaRp- 2 sin 2 9). (44) 
From Il2.'jll we then find the only non-vanishing components of the extrinsic curvature (up to symmetry) to be 

aM{3R 2 p 2 + a 2 [p 2 - 2 {R 2 + a 2 ) cos 2 9} } sin 2 6 

K m p 3 A i/2 S i/ 2 ( 45 ) 



and 



Ma 3 i?A 1 / 2 sin 3 cos6> , . 

Ke, = ~2 _5 . (46) 
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The electric part of the Weyl tensor now can be computed from © 

MR(2E + 3a 2 A sin 2 9)(p 2 - 4a 2 cos 2 6) 

>a 2 M(R 2 + a 2 )(3p 2 - 4a 2 cos 2 6) sin (9 cos6> 

MR(p 2 - 4a 2 cos 2 0)(£ + 3a 2 A sin 2 (9) 
7^ 



RR = (4,) 

= 3 ^ (48) 

£»e = ^ (49) 



MRY.ip 2 - 4a 2 cos 2 6) sin 2 6> 
E 4>4> = "g ( 50 ) 



and the magnetic part from Q 



aM(3p 2 - 4a 2 cos 2 60(21] + 3a 2 A sin 2 6) 
B RR = c °sv (5!) 



, aMR(R 2 + a 2 )(p 2 - 4a 2 cos 2 ■ 
oM(3p 2 - 4a 2 cos 2 0)(E + 3a 2 A sin 2 ■ 



-Brs = -3 — ; smO (52) 



Bee = cos6» (53) 

P £ 

aM£(3p 2 -4a 2 cos 2 0) . a 
5^ = v ^ 8 sur 6 cosO. (54) 

From and -By we construct Cy as well as the invariants / and J, following Section [HI As expected we find 

5=1, (55) 

which is by no means surprising but reassuring. 

B. Bowen— York Initial Data 

Initial Data describing rotating black holes can also be constructed by solving the constraint equations in the 
Bowen- York formalism. This approach assumes that the spatial metric is conformally flat. For rotating Kerr 
black holes, slices of constant Boyer-Linquist time are not conformally flat, nor are axisymmetric foliations that 
smoothly reduce to slices of constant Schwarzschild time in the Schwarzschild limit 01 • This suggests that Bowen- 
York initial data give rise to a spacctime that is distinct from the Kerr-Schild spacetime. This was demonstrated 
explicitly by Gleiser, Nicasio, Price and Pullin (|18|. hereafter GNPP), who used a perturbative calculation to show 
that Bowen- York initial data evolve into a spacetime that can be interpreted as a Kerr black hole with gravitational 
radiation. GNPP found that the amplitude of the emitted gravitational radiation scales with the square of the black 
hole's angular momentum L, and the power accordingly with L A . As an important test we verify in this Section that 
the BB scalar picks up this gravitational radiation and identifies the correct scaling. 

In a conformal transverse-traceless decomposition of the vacuum constraint equations the Hamiltonian constraint 
reduces to 



2 
3 

and the momentum constraint is 



8V 2 ip-ipR- -il; b K 2 + i>- 7 Ai^A 13 =0 (56) 



VjA* - \i) e f j VjK = (57) 

(see, e.g., 0,0^3 f° r recent reviews). The conformal factor tp relates the physical spatial metric to the conformally 
related metric 7y via 

lij =^%j, (58) 

where Vj and R are the covariant derivative and scalar curvature associated with jij , and the extrinsic curvature -Ky- 
is decomposed into its trace K and a traceless part Aij as 

Kij = ^ 2 A l} + pyK. (59) 
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Assuming conformal flatness (7^ = rjij , where rjij is the flat spatial metric in any coordinate system) and maximal 
slicing results in R = and K = 0, so that equations l|56|) and (|57|l reduce to 



1 A - A tJ 



(60) 



and 



(61) 



Here Vj is now the flat covariant derivative associated with rjij (as we assumed 7^ = rjij). 

The momentum constraint (|61|l decouples from the Hamiltonian constraint l|60() under these assumptions and can 
be solved analytically. For a spinning, unboosted black hole in polar coordinates (i.e. jij — diag(l, r 2 , r 2 sin 2 9)) the 
only non-vanishing component of Aij is 



.4 



3 sin 2 9 



r4> 



(62) 



where L is the angular momentum [l6| . 

Given the angular momentum L and Eq. 162|) , the conformal factor ip can be found from Eq. 1|60|) . Since the angular 
momentum L only enters squared in the Hamiltonian constraint l|6U|l . the conformal factor ip is even in L. In general, 
a solution to this quasi-linear elliptic equation has to be constructed numerically, but an approximate solution up to 
0(L 4 ) is given by 



M 
27 



where 



ipo 



-Y 

M 2 J 



M 
5r 



M 
27 



^0 + ^2 (cos 6>) +F{r,9)L 



, M\ 3 {M\ A (M 



(63) 



(64) 



and 



<T2 



1 (M 

10 \V 



(65) 



and where the P2(cos9) — (3cos 2 # — l)/2 is the second Legendre polynomial (see GNPP.) We allow for an un- 
known 0{L 4 ) contribution by including the unknown function F(r, 9) that, in principle, can be found by solving the 
Hamiltonian constraint to 0(L A ) in a method similar to that of GNPP. 

The extrinsic curvature, which must be odd in L, can then be found from Eq. (|59|) 



K r< p = ip A rcf> 
= 12 



sin 2 9 



(2r + M) 5 



48 [8(1 - 3 cos 2 9)r 3 
T 



20Mr 2 - 8M 2 r - M 3 ] sin 2 9 T , 



(2r + Mf 



L 3 + 0{L 5 ) 



(66) 



Even though we will find the BB radiation scalar £ is of order 0{L ), it turns out that it depends only on terms in ip 
up to order 0(L 2 ), so that it is independent of the unknown function F(r,9). A heuristic argument for this behavior 
can be given as follows (for a more rigorous proof see Appendix [DJ) The Bowen-York solution is equivalent to the 
Kerr solution up to O(L), and deviations enter only at order 0(L 2 ). In a transverse frame we therefore expect the 
leading-order terms in the Weyl scalars ipo and ip4, which are measures of the gravitational radiation that is absent 
in the Kerr solution, to be 0(L 2 ). The leading order term in the radiation scalar £, being the product of ipo and ipi) 
should therefore be 0(L 4 ), and higher order corrections in the conformal factor should enter at higher order. 

To construct £, we first compute the electric (01 and magnetic (0 parts of the Weyl tensor. We list the results in 
Appendix [CI and point out that these quantities do depend on F(r,9). We then compute the matrix C*- and, using 
© and ijIDfr. the invariants / <|G14() and J (|CJ 15|) . Again, I,, J depend on F(r,9). From ijTTJl we find the speciality 
index S to 0(i 4 ) 



S = 1 



33 2 10 ( 4r 2 _ 21Mj . + M 2)2 r 4 gin 4 

~~ 5 1 M 4 (2r + M) 12 



(67) 
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FIG. 2: The dimensionless BB scalar M 12 £/(L 4 sin 8) for Bowen-York initial data as a function of r/M. (For 6 — 0, 7r, the 
BB scalar £ = 0.) 



which turns out to be independent of F(r, 9). Because we were able to express the speciality index as a deviation from 
unity — i.e., we expressed spacetime as a deviation from algebraic speciality — we can find the BB scalar £ immediately 
from JI2J to 0(L 4 ) 



L (68) 



3 2 2 22 (4r 2 - 21Mr + M 2 ) 2 r w sin 4 6 
^~ 5 2 M 2 (2r + M) 24 

(see Fig.HJ). To 0(L 2 ) the Coulomb scalar X is given by 

Mr 3 r 4 768 

X = - 64 — r - 768 i cos 0- L + (24r 4 + 8r 4 cos 2 + 24r 3 M cos 2 

(2r + M) 6 (2r + M) 8 5 

- 8Mr 3 + 2r 2 M 2 cos 2 9 - 14M 2 r 2 - 8M 3 r - Af 4 ) — — ^ — L 2 . (69) 

; M(2r + M) 12 

As expected, the BB radiation scalar scales with L 4 . Since the deviation between Bowen-York data and a Kerr 
spacetime enters at order L 2 , the amplitude of the emitted gravitational radiation also scales with L 2 , and the power 
carried by the gravitational waves accordingly with L 4 (see GNPP.) As we have seen for the linearized waves on flat 
backgrounds in Section IlIII the BB radiation scalar scales with the square of the gravitational wave amplitude. At 
great distances (r 3> M), the BB scalar £ ~ L 4 /(M 2 r 10 ). The peeling off of the BB scalar with r~ 10 is understood 
as follows: The BB scalar is expected to drop like the product of V'o an d tp4- The usual peeling off of the product 
is expected to scale like M 2 /r 6 . However, in the Kerr spacetime at that order it would be multiplied by a zero 
coefficient. The next order, which would peel off like L 2 /r a , also vanishes because the deviation of the Bowen-York 
initial data from Kerr is at second-order in L in the metric functions. The next term is indeed the one found. This 
scaling is quadratic in x — XKcrr, as expected. When £ is calculated by finding the difference between / and x 2 5 the 
Kerr contributions to x cancel with the Kerr contributions to /, so that the BB scalar £ depends only on the radiative 
Bowen-York degrees of freedom. 

Notice, that £ peaks at r = M/2, i.e., at the Schwarzschild black hole's event horizon. The radiation described by 
the BB scalar inside the event horizon evidently cannot change the exterior spacetime. However, because the greatest 
deviations of the Bowen-York initial data from Kerr are localized near the event horizon, it is expected that much of 
this spurious radiation be absorbed by the black hole. 
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V. KASNER SPACETIME 



The Kasner spacetime is a homogeneous vacuum spacetime. From its spacetime metric we identify a unity lapse, a 
vanishing shift vector, and the spatial metric 

da 2 = t 2pi dx 2 + t 2p2 dy 2 + t 2p3 dz 2 , (70) 

where the parameters pi satisfy p\ + p2 + pz = 1 = p\ + p\ + p\. It is convenient express P2 and p^ in terms of p%, 
which we will call p, so that 

p = Pl p 2) a = (l-p±e)/2, (71) 

where g = [(1 — p)(l + Sp)] 1 ^ 2 . The Kasner spacetime is sufficiently simple to construct the BB radiation scalar in 
closed form, and also to demonstrate the construction of the quasi-Kinnersley frame from spatial data as they are 
used in numerical relativity. 

We compute the extrinsic curvature from the spatial metric and find 

K ab = diag{-pt 2 v- 1 ,-(l-p+Q)t-i>+e/2,-(l-p-Q)t-i>-e/2}, (72) 

so that K = —1/t. For the Kasner spacetime the magnetic part of the Weyl tensor vanishes, so that C a b = E a \,. The 
only non-vanishing components are 

C\ = p{l-p)t~ 2 

C\ = (l-p + g)(l+ P -g)t- 2 /A (73) 
C\ = (l-p-g)(l+p+g)t- 2 /A 1 

From © and IjlOfl we now find the Weyl tensor invariants 

/ = (74) 
J = 2f6 , (75) 

and from (|llfl the speciality index 

27 

S = -(l-p)p 2 . (76) 

Since the matrix C a b is diagonal, its three eigenvalues and corresponding eigenvectors can be found quite easily. In 
order to find which one of these corresponds to the quasi-Kinnersley frame we parametrize p as p = — 1/3 + g where 
q does not need to be small. Inserting this into l|7t)|l we find 

27 54 , 27 o 
S-l = -- q+ - q 2 -- q \ (77) 

The case q = corresponds to an algebraically special spacetime (Petrov type D), and for small values of q we can 
therefore use the expansion lfT^|) to find 

For small values of q this coincides with A x /2, so that we can identify 

We can find the BB scalar £ either from the series expansion l|17|) . or, having found \ already, from l|15fl . which yields 

5/1 o _\2 I 1 ^ „ \ _ 1 „2 /■„ „ \2 
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Evidently we have £ ^ 0, even though no radiation is present in the Kasner spacetime. [There are three cases for which 
£ = 0: q = (spacetime is Petrov type D), and q — 1/3 or 4/3, which are flat spacetimes (Petrov type 0).] While 
this seems very counter-intuitive, we remind the reader that the BB radiation scalar is an unambiguous measure of 
gravitational radiation only in a radiation zone, which the homogeneous Kasner solution does not possess. The name 
"radiation scalar" is therefore may be misleading in this case. 

Our results for /, J, and \ for the Kasner spacetime agree with those obtained through a pure spacetime approach 
by Cherubini et al. 20], who also compute the Weyl scalars ipo and tp4 in a transverse frame, and their product agrees 
with our value of £, as expected. 

We now demonstrate how the quasi-Kinnersley frame can be constructed from spatial data. The eigenvector that 
corresponds to the eigenvalue X X1 and hence to the quasi-Kinnersley frame, is 

a a = sS a x , (81) 

where s — t~ Pl is a normalization constant. Evidently er° is real and has vanishing imaginary part. 

Following paper I, the two real projections onto the hypersurface of the two real null vectors of the spacetime tetrad 
are 

A a = t~ pi 8 a x v a = -t~ pi 5 a x , (82) 

so that, with the unit normal to the hypersurface f° = d a , the two real spacetime null vectors of the frame are, up 
to a spin-boost parameter c, given by 

t = ^(d t + r^d x ) (83) 

■{d t -t-^d x ). (84) 



V2 



Notice, that the two real projections of £ a ,n a on the hypersurface are anti-parallel. This happens because by con- 
struction the normal f a lies in the spacetime tangent 2-plane spanned by l a and n a . In this degenerate case, the 
complex basis vector m a can be found as follows (see paper I for details) . We first choose an arbitrary real unit vector 
f a in the spatial 2-plane orthogonal to A° = — v a . Then, 



mr = -= (f a + i e abc X b f c ) . (85) 
We find that e abc — — [a b c]/ y/det(g mn ) = ~[ a b c]/t, with [a b c] being the permutation symbol. We next choose 



f a = t- p2 5%, so that 



m a = —( t- p2 d y - i r P3 d z ) . (86) 

The vectors I|83|l . (l84ll . and l|86|l are the vectors that make the quasi-Kinnersley frame. Note that we are able to 
determine the quasi-Kinnersley frame only up to spin-boost with parameter c = \c\ exp(ii9). In some situations one 
can choose the spin-boost parameter based on the physical properties of spacetime, so that the quasi-Kinnersley tetrad 
can be found (see, e.g., Ref. HJ). However, this problem awaits further study. 



VI. SUMMARY AND DISCUSSION 



The BB radiation scalar £ was introduced as an invariant measure of gravitational radiation in regions of spacetime 
where such radiation is unambiguously defined. However, it is uniquely and smoothly extended throughout a generic 
spacetime. In this paper, we have computed £ for a variety of analytical spacetimes and initial data sets, adopting a 
formalism that relies only on spatial data as they are typically used in numerical relativity. These calculations have 
illustrated the procedure used to define £ explicitly. That procedure has been described previously 0, 0, ^| , but not 
implemented. Its actual implementation here, therefore, should help clarify the procedure in general. 

For those examples we have considered here which unambiguously contain gravitational radiation, namely linearized 
waves on flat backgrounds and Bowen-York initial data for rotating black holes, our radiation scalar scales with the 
square of parameters that govern the gravitational wave amplitude. This suggests that it does indeed provide a 
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reasonable measure of the gravitational wave content of these data. (The speciality index S, in contrast, is indepen- 
dent of the waves' amplitude when the Coulombic part of curvature is perturbative.) We also have examined the 
homogeneous, cosmological Kasner solutions. In this case we found a non-vanishing £ even though Kasner does not 
contain gravitational radiation. We explain this apparent contradiction by noting the absence of a radiation zone 
in the Kasner spacetime. Thus, although the name "radiation scalar" for £ is misleading and inappropriate in this 
case, £ can still be found uniquely. From the mathematical point of view, it remains an interesting gauge-invariant 
"observable" even when it is not a radiation scalar. 

This paper has illustrated the potential of £ as a tool for discriminating among initial data sets for numerical 
relativity. However, it has done so using analytic, rather than actual numerical, data. An actual application to 
such numerical data sets would engender a fresh set of difficulties, originating in the fundamental discreteness of the 
variables involved and the associated subtlety in obtaining the required smooth extension of the quasi-Kinnersley 
frame into the strong-field regions of the initial data set. A full discussion of these issues lies outside the scope of this 
paper, and will be addressed elsewhere. However, we must make some comment here. To illustrate the resolution we 
envision to these difficulties in actual numerical work, we examine the results of Figure in the Appendix. These 
figures show all three values of the Coulomb and BB scalars, appropriately scaled to give dimensionless quantities, for 
linearized Einstein-Rosen waves of flat spacetime. At large radii, one branch of the upper graph, the dotted one, is 
clearly larger than the other two. This is the principal branch of £, and is associated with the quasi-Kinnersley frame 
in this region. Moving inward, there are three critical points where we must be careful, at p/a = y/2, 1 and l/y/2. 

The outermost critical point poses no real problem since the root we had been tracking, £ , is visibly discontinuous 
there. If we were moving inward in a numerical spacetime, calculating eigenvalues of the tensor C l j with sufficiently 

small radial steps, there is no way we would be confused about which eigenvalue was appropriate just inside p/a = y/2. 
Inside this first critical point, we would naturally be tracking the dashed root, rather than the dotted. 

The next critical point, at p/a — 1 poses a greater difficulty. The discontinuity there occurs in the derivative of the 
Coulomb scalar, rather than in the scalar itself, and could therefore be harder to detect in numerical data. However, 
even in this case, there is a simple resolution to the problem. We would be able to recognize numerically that the 
eigenvalue we are interested in, associated with the dashed curve, is degenerating with another, associated with the 
solid curve. There are then at least two things we could do to make sure the Coulomb scalar remains smooth at 
this critical point. One would be to use higher-order approximants, such as splines based on x-values at several grid 
points in each direction around the critical point, to explore which branch we ought to pick to preserve smoothness. 
An even simpler approach would be to exploit the known analytic structure of the Riemann surface underlying the 
multi-valued complex function underlying \- That structure is described in detail in Papers I and II. Essentially, this 
approach would boil down to calculating how the speciality index S varies in a neighborhood of the critical point, 
and mapping spacetime in that neighborhood into the Riemann surface. The critical point itself occurs where S lies 
exactly on the branch line we have chosen for \. As S crosses that branch line, there is no ambiguity at this level in 
how the branch of x ought to change to preserve analyticity. Thus, once we have detected the possible presence of 
such a critical point in numerical data, we could use the explicit formula of 11211 and the local structure of the function 
S on spacetime to predict the exact eigenvalue which corresponds to the quasi-Kinnerlsey frame at every grid-point 
nearby. 

The third critical point, at p/a = l/y/2 poses no difficulty whatsoever. At this point, the quasi-Kinncrslcy value 
is associated with the solid curve, and it is the other two which are discontinuous. We would scarcely notice this in 
actual numerical work. 

The Einstein-Rosen example discussed here is actually quite useful. The three critical points we have observed 
in this example typify the only three kinds of critical points which can occur in a general spacetime, whether nu- 
merical or analytic in origin. The above discussion shows how these three types of critical point may be handled 
using a combination of simple numerical tests and insights from the analytical development of the radiation scalar 
approach. Although further testing will doubtless be required in developing general numerical implementations, these 
observations are encouraging for the ultimate viability of this approach in numerical relativity. 

A second question which lies outside the scope of the present paper, but which we nonetheless ought to address 
briefly here, concerns how the radiation scalar could potentially be used in discriminating among potential choices 
of initial data. Consider, for example, initial data describing a compact binary at small binary separation. As 
we have discussed before, these data will presumably contain "astrophysically sound" radiation together with the 
"junk radiation" that we would like to minimize. Such initial data can be constructed using different approaches, 
namely different decompositions of the constraint equations and different choices for the freely specifiable background 
geometry. Different choices lead to physically distinct solutions 0,0,0 with distinct gravitational wave content. Our 
comparison of Kerr and Bowen-York data for rotating black holes in Section IIVI illustrates how the BB radiation 
scalar may be used as a diagnostic for comparing the amount of gravitational waves that these initial data sets would 
emit when evolved dynamically. 

We find in the analytical examples considered here that the BB radiation scalar scales as expected with a parameter 
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that controls the gravitational wave content of these examples. Such a parameter may not exist in numerically 
generated initial data sets, in which case it may be less clear how any non-zero value of £ should be interpreted. 
The BB radiation scalar may nevertheless be a useful diagnostic. Even for a single initial data set the ratio £/x 2 
would provide a dimensionless measure of the strength of the "radiative" fields compared to the "coulombic" fields. 
Perhaps more importantly, it may be useful to compare £ for different initial data sets that approximate the same 
physical situation. For example, such a comparison may provide some guidance for deciding which choices in the 
construction of initial data sets yields a more truthful representation of binary black holes or neutron stars at small 
binary separations. We can envision various different ways of how such a comparison could be made. One potentially 
useful approach would focus on the radiation zone, where we have reliable intuition about how the real radiation 
content ought to fall off. Comparing the fall-off behavior with that of the BB radiation scalar may help distinguish 
this "astrophysically sound" radiation content from the "junk" content and minimize the latter. 

As a word of caution, however, we point out that such a comparison can at best provide some guidance rather 
than conclusive evidence. It is possible that the individual radiation fields ipo and ^4 change in such a way that their 
product £ = ipoipi decreases even though the true radiation content increases. However, for many generic sets (or 
families) of initial data this may not be the case. Since a diagnostic providing conclusive evidence is not yet available, 
we therefore believe that the BB radiation scalar may in the meantime provide some useful guidance as long as we 
keep these caveats in mind. 

We finally discuss some further limitations of this approach. While the BB radiation scalar £ contains some measure 
of gravitational wave content, it is not clear immediately how this information can be translated into gravitational 
wave templates that might be useful for gravitational wave observers. Specifically such observers would presumably 
measure the gravitational wave amplitudes h x and h + in a transverse-traceless frame. In the linearized examples 
we have studied here, both the intuitive gravitational wave content and the BB radiation scalar scales with the 
square of a parameter describing intuitively the strength of the waves being modeled. This does not mean, however, 
that £ is proportional to the energy of the gravitational wave. In general the latter scales with \i))q + 1/^4 1 2 in the 
appropriate basis, while £ = ipo4>4- Since ipo and -04 have different fall-off behavior, £ also falls off differently from 
both the gravitational wave amplitude or energy. In short, although £ describes only the radiative degrees of freedom, 
£ does not carry with it all the information — or even sufficient information — about the gravitational waves required 
by observers. This issue will be discussed more fully in a forthcoming paper. 
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APPENDIX A: DIRECT CALCULATION FOR EINSTEIN-ROSEN WAVES 

In this appendix we find x an d £ for Einstein-Rosen cylindrical waves directly from Eqs. I|12|) and (|13|) . First, we 
re-write Eq. l(T2|) as 



where 



Vs - vs^i 



1/3 



-1/3 



(Al) 



(A2) 



Here, a is one of the three cubic roots of unity, «o — U = exp (2iw/3), and a_ = exp (4«7r/3). By the cubic and 
square roots in Eq. I|A2(I we mean the principal branch roots for either. The choice of the cubic root then is done 
through the choice of on. This choice labels the three different branches of the function Z % {S). The corresponding 
functions £ 0,± can be calculated using 



e>± = / -3(x°- ± ) 2 



(A3) 



or, equivalently, from Eq. I|13fl using a similar approach to its three branches. 
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p/a 



FIG. 3: The three branches of the Coulomb scalar x '* (top panel) and the BB scalar £ 0,± (bottom panel) for linearized 



Einstein-Rosen cylindrical waves. The principal branch — \ an d £ 
the solid curves. 



-are dotted, \ an d ip are dashed, and \ an d ip are 



Figure shows the three different branches for x 0,± and £ 0,± , correspondingly. Notice that the three branches for 
either function are non-differentiable, and for x° and x~ even discontinuous. For either function, however, branches 
can be changed so that the resulting functions are smooth. The identity of the physically meaningful functions can 
be determined by demanding a proper asymptotic behavior of x an d £■ In this case, starting at great distances, the 
chosen branch is the principal branch, x° (and correspondingly, At the branch point p/a = y/2 we change the 
branch to x~ (and, correspondingly, and at p/a = 1 we change the branch to x + (and, correspondingly, £ + ). 
The resulting smooth curves are identical with those presented in Fig. ^ 

APPENDIX B: SOME EXPRESSIONS FOR LINEARIZED QUADRUPOLE WAVES 



From the spatial Ricci tensor we find that at the moment of time symmetry, t = 0, the only non- vanishing 
components of C l j are to leading order in e, 



C\ = 24(25r b A 3 -2r 8 A 4 + 54r^A-78r 4 A 2 + 224r 4 A 2 cos 2 0-73r b A 3 cos 



C", 



c\ = 



^158 r 2 A cos 2 9 + 6 r 8 A 4 cos 2 



4 cos 2 6) e- r A 



c/, 



~ 2 C r e = 24 sin 9 cos 6* (-12 + 293 r b A - 64 r 8 A + 150 r A + 4 r ' A - 420 r A ) e~ r A e/r 
-24 (6 + 715 r 6 A 3 - 290 r 8 A 4 + 78 r 2 A + 292 r 8 A 4 cos 2 9 - 739 r 6 A 3 cos 2 9 
-567 r 4 A 2 + 640 r 4 A 2 cos 2 9 + 43 r 10 A 5 - 43 cos 2 6 r w A 5 — 130 r 2 A cos 2 6 



-4 cos 2 9 ~ 2 r 12 A 6 + 2 cos 2 6 r 12 A 6 ) 



-r A 



e/r 2 



24 (6 + 690 r 6 A 3 - 288 r 8 A 4 + 24 r 2 A + 286 r 8 A 4 cos 2 6 - 666 r 6 A 3 cos 2 6 
-489 r 4 A 2 + 416 r 4 A 2 cos 2 9 + 43 r 10 A 5 - 43 cos 2 r 10 A 5 + 28 r 2 A cos 2 (9 



-8 cos 2 9 - 2 r 12 A 6 + 2 cos 2 r 12 A 6 ) 



-r A 



e/r\ 



(Bl) 
(B2) 

(B3) 
(B4) 
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Restricting analysis to the equatorial plane, 8 = tt/2, we find the curvature scalars 

/ = 576 (36 - 1548 r 4 A 2 - 51744 r 6 A 3 + 353559 r 8 A 4 - 27664 r 18 A 9 - 172 r 22 A 11 - 451732 r 14 A 7 + 612 r 2 A 

-773889 r w A 5 + 803755 r 12 A 6 + 3005 r 20 A 10 + 146051 r 16 A 8 + 4 r 24 A 12 ) e~ 2 r * A e 2 /r 4 (B5) 

and 

J = -6912 A (2 r 6 A 3 - 54 + 78 r 2 A - 25 r 4 A 2 ) 

x (-6 - 690 r 6 A 3 + 288 r 8 A 4 - 24 r 2 A - 43 r 10 A 5 + 489 r 4 A 2 + 2 r 12 A 6 ) 

x (-6 - 715 r 6 A 3 + 290 r 8 A 4 - 78 r 2 A - 43 r 10 A 5 + 567 r 4 A 2 + 2 r 12 A 6 ) e~ 3r2 A e 3 /r 4 (B6) 
From these we compute the speciality index 

S = 27 r 4 A 2 (2 r 6 A 3 - 54 + 78 r 2 A - 25 r 4 A 2 ) 2 (-6 - 690 r 6 A 3 + 288 r 8 A 4 - 24 r 2 A - 43 r 10 A 5 + 489 r 4 A 2 + 2 r 12 A 6 ) 2 
x (-6 - 715 r 6 A 3 + 290 r 8 A 4 - 78 r 2 A - 43 r 10 A 5 + 567 r 4 A 2 + 2 r 12 A 6 ) 2 

/ [4(36 - 1548 r 4 A 2 - 51744 r 6 A 3 + 353559 r 8 A 4 - 27664 r 18 A 9 - 172 r 22 A 11 - 451732 r 14 A 7 + 612 r 2 A 

- 773889 r 10 A 5 + 803755 r 12 A 6 + 3005 r 20 A 10 + 146051 r 16 A 8 + 4 r 24 A 12 ) 3 ] . (B7) 

As expected, S is independent of the wave amplitude e. Note, that for Ar 2 -» oo, 5^0. This is in spite of the 
fact that spacetime approaches algebraic speciality. In fact, spacetime approaches algebraic speciality exponentially 
fast, as both / and J decay exponentially with Ar 2 . This behavior is captured by the BB scalar £: from Eq. Ijl3|l . £ 
equals / times a function of S only. The speciality index S is dropping off at great distances to zero like A -4 r -8 , so 
that the function that multiplies / approaches a constant. The BB scalar £ drops off exponentially, specifically like 
exp(-2Ar 2 ). 

Figure 2] shows the speciality index S as a function of X 1 / 2 r on the equatorial plane. There are 15 values of r for 
which 5=1, arranged in five triplets. The values of p = A : / 2 r for which S — 1 arc 

iPi,2,3 = 0.30800524515557994139, 0.45049403784384619791, 0.59829320761717622144 

2 pi! 2 . 3 = 1.0335527869353531668, 1.0467779593567221747, 1.0682855799176748686 

spi^a = 1.6984496413703925629, 1.7168422939913076313, 1.7311171061162415883 (B8) 

4Pi,2.3 = 2.4478787994829345025, 2.4658595692782475639, 2.4824885271217855892 

5Pi^3 = 3.3409626177591728790, 3.3526507427990501828, 3.3643720160172523032 

In between points at which S = 1, we have points at which \S — 1| = 1, which occur at the following locations: 

lP * 1 23 = 0, 0.39945741403776425339, 0.50160586730815486639 

2Pi'2' 3 = 0.98542263099668675629, 1.0416944916187449946, 1.0527498193777135990 

3^ 3 = 1.7112566884721670706, 1.7219786282054512104 (B9) 

4P123 = 1.8534934086232694765, 2.4600232831874878020, 2.4715457963797601620 

5 pl^ 3 = 2.8449085086360266876, 3.3487504415432614190, 3.3565547217000555743 

Notice that there is no value for 3/?*: the speciality index S does not vanish between 2P3 and 3 p\. In addition, 
\S - 1| -> 1 also as A 1/2 r -> 00. 

It appears that in some cases - specifically for odd-parity quadrupole waves - the curvature invariant J may vanish 
to 0(e 3 ), so that its leading order is 0(e 4 ). Generic gravitational waves include both polarization states, so that the 
total curvature invariant J is still of 0(e 3 ), and S — O(e ). 

APPENDIX C: RESULTS FOR BOWEN-YORK INITIAL DATA 

Inserting Q66J) into 

Bij = -e ( fcZ V fc X y = -7 im e mkl (K ljtk - K ln T] k ) (CI) 
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FIG. 4: The speciality index S as a function of A 1 / 2 r on the equatorial plane, for linearized, even-parity, quadrupole waves. 
The top panel magnifies the regions where S oscillates rapidly. 



and symmetrizing we find that the only non- vanishing coefficients of the magnetic part of the Weyl tensor are 

64 cos 9 r 



(2r + M) 4 



1536 [16(2 - 3 cos 2 9)r 3 - 2QMr 2 - 8M 2 r - M 3 ] cosfl 3 
+ "To M(2r + M) 10 ( ' 

48 r (2r - M) sin 9 
Bre (2r + M)5 ^ 

384 [40(1 - 3 cos 2 6»)r 4 - 20(7 - 3 cos 2 9)Mr 3 - 16M 2 r 2 + 4M 3 r + M 4 ] r sin 9 3 
+ 1 (2r + M)" ^ (C3) 

Bfle = 48 COS - £ 
0e (2r + M) 4 

384 cos<9 [20 Mr 2 + 8 M 2 r + M 3 + 24 r 3 cos 2 9- 8 r 3 ] r 2 3 

+ ~5 (2r + M) 10 M ( } 

cos 6* sin 2 6* r 2 

^ = 48 (2r + M) 4 L 



384 cos (9 sin 2 (9 [72 r 3 cos 2 (9 - 56 r 3 + 20 Mr 2 + 8 M 2 r + M 3 ] r 2 . 
+ 1 (2r + M)i"M ^ (L5) 

The electric part of the Weyl tensor is computed from the Ricci tensor, whose non-vanishing components to 0(L 4 ) 
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arc 



8M 
(2r + M) 2 r 



32 144 r 4 cos 2 9 — 360 r 3 M cos 2 + 36 r 2 M 2 cos 2 — 32 M 2 r 2 - 48 r 4 — 8 M 3 r + 440 Mr 3 — M 4 
+ T (2r + M) 8 Mr J 

64 

+ — [-8896 Mr e + 16176M 2 r 5 + 6480 M 3 r 4 + 800 M 4 r 3 - 36 M 5 r 2 - 12 r M 6 - M 7 + 960 r 7 
25 

- 10704 r 5 cos 2 M 2 + 31872 r 6 cos 2 M - 14400 M r 6 cos 4 + 2016 M 2 r 5 cos 4 - 240 M 4 r 3 cos 2 6 
+ 36 M 5 r 2 cos 2 - 3360 r 4 cos 2 M 3 - 5184 r 7 cos 2 + 8064 r 7 cos 4 0] L 4 A(2 r + M) 14 M 2 r] 



4L 4 |"2M(2r-M) cot0 Q 2 2r + 3M 1 
2r + M [ r(2r + M) 2 r~ 8 ~ 2r + M r r 



- 2rd r 



F(r, 



(C6) 



r 2 (2r-M) 2 

— sin cos 0— jX _ — L 

5 (2r + M) 7 M 



+ 



768 
~T 
1536 
"25" 



r 2 sin 9 cos 6» (156 r 4 cos 2 6» - 52 r 4 - 78 r 3 M cos 2 6» + 186 M r 3 + 26 M 2 r 2 



3 M r — M ) L / [(2 r + M) M ] + 



2r + M 



L 



M 



2r + M 



de- rdre F(r,9) 



(C7) 



Ree 



4Mr 



(2r + M) 2 

16 (168 r 4 cos 2 9 - 72 r 4 - 304 Mr 3 + 384 r 3 M cos 2 + 42 r 2 M 2 cos 2 - 38 M 2 r 2 -8M 3 r-M 4 )r 2 



(2r + M) 8 M 

- 42AfV - 12M 6 r + 8M 4 r 3 + 13284r 7 cos 4 

- 11136 r 7 cos 2 + 8640 r 6 M cos 4 + 3360 r 6 M cos 2 + 3456 r 5 M 2 cos 4 + 2880 r 5 M 2 cos 2 + 2736 r 4 M 3 cos 2 
+ 552 r 3 M 4 cos 2 + 42 r 2 M 5 cos 2 9)r L 4 / [M 2 (2r + M) 14 ] 



32 

— (1152r 7 - M 7 - 3424Mr 6 + 1152M 2 r 5 + 384M 
25 



4r 



2r + M 



L 



M(2r-M) „ r(6r + M) „ „ , „ 

j + COt09 e + ^— — -i- 8 r + 2dgg + T 2 8 r 



(2r + M) 2 



2r + M 



F(r,9) 



(C8) 



Rd 



AMr 



sin 



(2r + M) 2 
16 

y 

32 
25 



sin 



2 (336 r 3 M cos 2 + 30 r 2 M 2 cos 2 + 120 r 4 cos 2 - 24 r 4 - 8 M 3 r - 256 M r 3 - 26 M 2 r 2 - M 4 ) r 2 



M{2r + Mf 

sin 2 0(1440 M 2 r 5 cos 4 + 576 Mr 6 cos 4 + 5760 r 7 cos 4 + 408 M 4 r 3 cos 2 9 

- 2688 r 7 cos 2 + 30 M 5 r 2 cos 2 + 3648 r 5 cos 2 M 2 + 2064 r 4 cos 2 M 3 + 10848 r 6 cos 2 M + 1056 M 3 r 4 
+ 152 M 4 r 3 - 2848 M r 6 + 2400 M 2 r 5 + 768 r 7 - 30 M 5 r 2 — M 7 — 12 rM 6 )rL 4 I [M 2 (2 r + M ) 14 ] 



4r 



2r + M 



sin 2 0L 4 



M(2r - M) 
(2r + M) 2 



+ r 2 9 rr + 2 cot <9 e + 



r(6r + M) 
2r + M 



d r + 9ge 



F(r,9), 



(C9) 
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and the nonzero components of the electric part of the Weyl tensor are given by 
8M 



Err ~r (2r + M) 2 



+ y 



32 144r 4 cos 2 + 36r 2 M 2 cos 2 - 48r 4 - 32 M 2 r 2 + 80Mr 3 - 8 M 3 r - M 4 



Mr(2r + M) 8 

64 

— (704 M r 6 - 7824 M 2 r 5 — 3120 M 3 r 4 - 400 M 4 r 3 - 36 M 5 r 2 - 12 r M 6 - M 7 + 960 r 7 
25 

+ 9792 r 7 cos 4 - 6912 r 7 cos 2 - 8256 r e cos 2 M + 16128 M r 6 cos 4 

+ 2448 M 2 r 5 cos 4 + 36 M 5 r 2 cos 2 + 960 M 4 r 3 cos 2 + 12864 r 5 cos 2 9 M 2 

+ 6240 r 4 cos 2 9 M 3 ) L 4 j[M 2 r{2r + M) 14 ] 

4 L 4 

+ 3 2r + M [ r(2r + M) 2 



6M(2r-M) n n 1 2r-5M„ cot 9 

i - 2r 9 rr + - dee + 2 d r + 

; - 2r + M r 



F(r,6) 



(C10) 



768 r 2 (2r-M) 

^re = -=- sin cos — tttttt 

5 (2r + M) 7 M 



+ 



1536 
~25~ 



sin 6» cos 6» r 2 (156 r 4 cos 2 - 52 r 4 - 78 r 3 M cos 2 



+ 186 Mr 3 + 26 M 2 r 2 - 3 M 3 r - M 4 ) L 4 / [(2 r + M) 13 M 2 } 



+ 4- 



M 



2r + M V 2r + M 



rc> re F(r,9) 



(Cll) 



4Mr 



(2r + M) 2 

16 (168 r 4 cos 2 - 72 r 4 - 304 M r 3 + 384 r 3 M cos 2 + 42 r 2 M 2 cos 2 9 - 38 M 2 r 2 - 8 M 3 r - M 4 )r 2 



5 

32 
25 



M{2r + Mf 

(23456 Mr 6 - 66048 M 2 r 5 - 26496 M 3 r 4 - 3352 M 4 r 3 - 42 M 5 r 2 - 12 r M 6 - M 7 

+ 1152r 7 + 13824r 7 cos 4 9 - 11136r 7 cos 2 9 + 29616 r 4 cos 2 9 M 3 
+ 89280 M r 6 cos 4 + 3456 M 2 r 5 cos 4 + 42 M 5 r 2 cos 2 

- 104160 r 6 cos 2 9M + 3912 M 4 r 3 cos 2 (9 + 70080 r 5 cos 2 9 M 2 )r L 4 l[M 2 {2 r + M ) 14 ] 



4 

+ - 



3 2r + M 



L 



cot 9 de 



r(2r - 5M) 
2r + M 



8 r -2dee+ r z 0, 



2... 3 M(2r-M) 



(2r + Mf 



F(r, 9) 



(C12) 



4Mr 2 
^ - (2r + M) 2 Sm ^ 

- ^ sin 2 0r(-384 r 3 M cos 2 9 + 30 r 2 M 2 cos 2 

5 

+ 120 r 4 cos 2 - 26 M 2 r 2 - 24 r 4 - 8 M 3 r + 464 Mr 3 - M 4 ) L 2 / [M (2 r + M) 8 ] 
32 

- — sin 2 r(-57024 M r e cos 4 + 5760 r 7 cos 4 + 1440 M 2 r 5 cos 4 
25 

- 2688 r 7 cos 2 - 17136 r 4 cos 2 M 3 - 44352 r 5 cos 2 M 2 + 30 M 5 r 2 cos 2 

- 1992 M 4 r 3 cos 2 + 87648 r 6 cos 2 M + 20256 M 3 r 4 - 22048 M r 6 + 768 r 7 

- 30 M 5 r 2 - M 7 + 2552 M 4 r 3 + 50400 M 2 r 5 - 12 r M 6 ) i 4 / [M 2 (2 r + M) 14 ] 



3 2r + M 



sin 2 0L 4 



3M(2r - M) 
(2r + M) 2 



dee+ r( \ ■ d r + 2 C ol9d e 
2r + M 



F(r,9). (C13) 
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Next, we find the tensor C a b to 0(L 4 ), from which we construct the Weyl tensor invariants I and J. To 0(L 4 ), we 



have 



12288 M 2 r b 
(2r + M) 12 



294912 i 



cos 9Mr 7 
(2r + M) 14 



294912 (128 r 4 cos 2 6> + 24 r 4 + 144 r 3 M cos 2 6> - 8 Af r 3 + 32 r 2 Af 2 cos 2 8 - U M 2 r 2 - 8 M 3 r - M 4 )r 6 



5 

1179648 



(2r + M) 18 

i cos 9 (80 r 4 cos 2 6>- 128 r 3 M cos 2 (9 + 20 r 2 M 2 cos 2 6> + 16 r 4 + 112 Mr 3 



136 M 2 r 2 -7M 4 



56 Af 3 r) r 7 L 3 /[(2 r + Af ) 20 M] 
196608 



25 



(300 M 7 r + 19 Af 8 - 32288 M 2 r 6 + 29712 M 3 r 5 + 21312 Mr 7 + 12972 Af 4 r 4 



+ 1548 Af 6 r 2 + 4672 M 5 r 3 + 2496 r 8 - 1920 r 8 cos 2 9 + 8640 r 8 cos 4 

+ 43776 r 7 M cos 2 - 98328 r 4 M 4 cos 2 (9 - 158016 r 7 M cos 4 9 

- 39504 r 5 Af 3 cos 4 9 - 1356 r 2 Af 6 cos 2 (9 + 540 r 4 M 4 cos 4 

- 196608 r 5 M 3 cos 2 6 - 33408 r 6 Af 2 cos 2 8 - 233184 r 6 M 2 cos 4 (9 

- 19056 r 3 Af 5 cos 2 (9) r 6 L 4 /[(2 r + Af) 24 AT 2 ] 

Mr 6 T/1 r . „„ „ r(2r - 5M) 



4096 



(2r + M) n 



L 4 



cot So 



2r + M 
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APPENDIX D: THE BOWEN-YORK BB RADIATION SCALAR IS INDEPENDENT OF F(r,0) 

In this Appendix we outline a calculation that shows why the speciality index S, and hence the BB radiation scalar 
£°, is independent of the 0(L 4 ) terms in the conformal factor. 
Let us expand the curvature invariants I and J to 0(L 4 ) as 

4 4 

f = ]T ifc L fc , J = Y / JkL k , (Dl) 

k=0 k=0 

where the expansion coefficients ik and jk are given by l|C14(l and l|C15|l . The only coefficients that depend on F(r, 9) 
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are in and j'4. The speciality index S is given by (|llf> and can also be expanded in powers of L up to order 0(L A ) as 

4 

S = ^Ts fcj L fc . (D2) 

fe=0 

The expansion coefficients can be found from the definition of 5* and the expansion coefficients for / and J. 

The leading order term in S is unity. A straightforward calculation shows that the next contributing order is at 
0(L 4 ). What remains to be shown is that the 0{L A ) term S4L 4 is independent of F(r, 8). It turns out that the only 
dependence of S4 on either 14 or j'4 is through the combination 2i j4 — 3j'o*4- A direct substitution finds that this 
combination equals exactly zero, which completes the proof. 

The BB scalar £° is proportional to / times a power series in S — 1, the leading order term being of 0(S — 1). 
Consequently, the earliest order at which F(r,9) can contribute is at 0(L 8 ). Specifically, it does not contribute at 
0(L 4 ), which is the leading order for the BB scalar. 
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